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Abstract 



> 
(N 

' Weak measurements are a new tool for characterizing post-selected quantum sys- 

tems during their evolution. Weak measurement was originally formulated in terms 
\ of von Neumann interactions which are practically available for only the simplest 

|0 I single-particle observables. In the present work, we extend and greatly simplify a 

recent, experimentally feasible, reformulation of weak measurement for multipar- 



ticle observables [Resch and Steinberg (2004, Phys. Rev. Lett., 92, 130402)]. We 
also show that the resulting "joint weak values" take on a particularly elegant form 
P5 ■ when expressed in terms of annihilation operators. 

O^' Key words: Weak measurement. Entanglement, Post-selection, Annihilation 

Operator 
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Weak measurement was originally proposed by Aharanov, Albert and Vaidman 
(AAV) [1] as an extension to the standard von Neumann ("strong") model of 
measurement. A weak measurement can be performed by sufficiently reduc- 
ing the coupling between the measuring device and the measured system. In 
this case, the pointer of the measuring device begins in a state with enough po- 
sition uncertainty that any shift induced by the weak coupling is insufficient to 
distinguish between the eigenvalues of the observable in a single trial. While 
at first glance it may seem strange to desire a measurement technique that 
gives less information than the standard one, recall that the entanglement gen- 
erated between the quantum system and measurement pointer is responsible 
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for collapse of the wavefunction. Furthermore, if multiple trials are performed 
on an identically-prepared ensemble of systems, one can measure the average 
shift of the pointer to any precision - this average shift is called the weak value. 
A surprising characteristic of weak values is that they need not lie within the 
eigenvalue spectrum of the observable and can even be complex [2,3,4]. On 
the other hand, an advantage of weak measurements is that they do not dis- 
turb the measured system nor any other simultaneous weak measurements or 
subsequent strong measurements, even in the case of non-commuting observ- 
ables. This makes weak measurements ideal for examining the properties and 
evolution of systems before post-selection and might enable the study of new 
types of observables. Weak measurements have been used to simplify the cal- 
culation of optical networks in the presence of polarization-mode dispersion 
[5], applied to slow- and fast-light effects in birefringent photonic crystals [6], 
and bring a new, unifying perspective to the tunneling-time controversy [7,8]. 
Hardy's Paradox, introduced in Ref. [9], was analyzed in terms of weak values 
in [10]. In Ref. [11], weak values were used to physically explain the results of 
the cavity QED experiment described in [12]. The opposing views expressed 
in Refs. [13,14] on the role of which-path information and the Heisenberg un- 
certainty principle in the double-slit experiment are reconciled with the use of 
weak values in [15]. Weak measurement can be considered the best estimate 
of an observable in a pre and post-selected system [16]. 

The von Neumann interaction was originally used to model standard quan- 
tum measurement by mathematically describing the coupling between the 
measured system and the measurement pointer [17]. The interaction couples 
an observable A of the quantum system to the momentum P of the pointer, 

n = gAP, (1) 

where g is the coupling constant which is assumed to be real to keep Ti Her- 
mitian. Since A and P act in different Hilbert spaces we can safely assume 
they commute. This interaction would be difficult to implement were it not 
for the fact that typically the measured system itself is used as part of the 
measurement device. When measuring A oi a particle an independent degree 
of freedom of the particle can be used as the pointer. For example, a bire- 
fringent crystal can be oriented so that it will displace the position of photon 
by an amount that depends on the photon's polarization [18]. Here, A is the 
polarization observable and the pointer is the position of the photon. Another 
example is the Stern-Gerlach apparatus, where A is the spin of the particle and 
the pointer is the momentum of the particle. If such a measurement strategy 
were not availabe, one would require a strong controllable interaction between 
the quantum system and a separate pointer system. This is typically far too 
technically difficult to implement. 

In modern quantum mechanics, we are increasingly interested in a different 
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class of observables than in the above example, in which only a single particle 
is involved. Often, one would like to measure correlations between observables 
of distinct particles, like S1S2, the spin of particle one times the spin of particle 
two. Moreover, any experiment that utilizes or directly measures properties 
of entanglement is based on such observables and so, much of quantum infor- 
mation and quantum optics deal with these composite or joint observables. 
The exciting results and complex, rich range of features discovered by studies 
of entanglement suggests that weak measurement of joint observables should 
also produce valuable and interesting results. In fact, there already exist a few 
theoretical ideas for weak measurements that center around joint observables, 
such as Hardy's Paradox [10], nonlocality of a single particle [8], and exten- 
sions of the Quantum Box Problem [19,20]. We call the weak value of a joint 
observable the "joint weak value." If the composite observable is a product 
of N single particle observables then the weak value is called the 'Wth-order 
joint weak value" . 



Joint observables are extremely difficult to measure directly with either strong 
or weak types of measurement. The difficulty lies in the fact that the nec- 
essary von Neumann interaction couples two separate observables, and hence 
particles, to a single pointer. One, therefore, can no longer use the extra de- 
gree of freedom on one of the particles as the pointer and so, one requires 
multiparticle interactions. An approach using multiparticle interactions was 
outlined in a proposal for a weak measurement experiment with ions but so 
far there have been no experimental weak measurements of joint observables 
[21]. On the other hand, experimental strong measurements of joint observ- 
ables are feasible and even commonplace. This is made possible by employing 
a different measurement strategy. Instead of measuring the joint observable 
directly, each single particle observable is measured simultaneously but sepa- 
rately. For example, instead of measuring S1S2 directly we can measure Si 
and ^2 separately and then multiply the results trial by trial. If one wants to 
strongly measure the joint observable A1A2...AN — M, instead of using the 
multiparticle von Neumann Hamiltonian Ti, = gMP, the general strategy is 
to simultaneously apply N standard single-particle von Neumann interaction 
Hamiltonians, 



H = giAiPi + g2A2p2 + ... (2) 

N 

-Y.9jAjPy (3) 

Given that we can already perform each of the single-particle Hamiltonians, it 
is straightforward to implement the total Hamiltonian. This strategy allows 
one to make projective measurements of M which is all that is required to 
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measure the expectation value of M, 

(m) = (iii2...iiv) oc {x^X2...Xr,) , (4) 

where Xi is the position operator of the pointer and provided all Ai commute. 
In other words, the expectation value of M is related to the correlation between 
the positions of all N pointers. 

In two earlier works, an analogous strategy was applied to weak measurements 
[22,23]. The Hamiltonian in Eq. (3) is utilized in the weak regime to create 
correlations in the deflections of the pointers proportional to the weak value. 
Specifically, the A^th-order joint weak value was related to two correlations 
between all N pointer deflections and a complicated combination of lower- 
order joint weak values. In this work, we show that the A^th-order joint 
weak value takes on an elegant and simple form closely related to the strong 
measurement formula in Eq. (4) when expressed entirely in terms of A^-pointer 
correlations. This new and simplified form lends itself to a new way of thinking 
about single and joint weak measurements in terms of expectation values of 
products of annihilation operators. 

We begin by deriving AAV's formula for the weak value of a single particle 
observable. AAV based weak measurement on the weak limit of the standard 
approach to measurement. Specifically, they use the von Neumann interaction 
in Eq. (1), which we assume to be constant over some interaction time t. The 
measurement pointer is initially in a Gaussian wavefunction centered at zero, 

where a is the rms width of |0(a;)|^. In most experiments, quantum me- 
chanical systems are initially prepared in a known initial state |/). Since this 
preparation usually involves measuring an ensemble of systems and selecting 
the subensemble with the correct outcome, this is called pre-selection. For a 
strong measurement, the von Neumann interaction with a pre-selected system 
state shifts the mean position of the pointer <^X^ by gt {I\ A \I) and leaves 

^P^ unchanged. AAV considered the case where we further restrict ourselves 
to the subensemble of system states that are found to be in \F) after the mea- 
surement, a procedure called post-selection. A weak measurement performed 
between the pre and post-selection can result in very different expectation 
values than in strong measurements, as we will see. 

After the pointer weakly interacts with the initial sytem-pointer state |'?/'(0)) = 
|/) the state evolves to 
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We project out the part of the state that is post-selected in state \F) , 
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(6) 
(7) 

(8) 



This leaves the state of pointer after the interaction and post-selection. In 

the limit of an ideal weak measurement, gt 0, |(-F|/)|^ =Probsucccss is the 
probability the post-selection succeeds [7]. If we renormalize the state and 
then truncate the amplitude of each term to lowest order in gt we get 



igt (F|i|J) ^, 



(9) 



which is just equivalent to dividing by {F\I) = i/Probguccess- The subscript 
fi, corresponding to final state \F) and intial state |/), labels the final pointer 
state, with which we can now calculate the expectation value of X of the 
pointer. The terms which contain an expectation value of an odd number of 
operators go to zero since the pointer is initially an even function about zero. 
To first order in gt, the remaining terms give us 
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Here, Qj- is used to signify the expectation value of a pointer observable only 
in the subensemble of measured systems that start in state |/) and are later 
post-selected in the state \F) . Similarly, the momentum expectation value is 
given by 
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The shifts from zero in both the X and P expectation values are proportional 

!Ctively, c 

{F\A\I) 



to the real and imaginary parts, respectively, of the weak value \^/^ which 
is defined as 



In fact, AAV showed that for sufficiently weak coupling the final 

pointer state, will be (^/2iTa^ ^ exp ^— (^x — (^A^^^ /4cr^^, unchanged ex- 
cept for a shift by the weak value. 

It has been argued that it is the backaction of the measurement on the mea- 
sured system that leads to a finite lm(^Aj^ and thus a nonzero (-P)g [7]. 

In addition, as the measurement becomes weaker (-P)^. becomes more and 



more difficult to determine; (-P/^. decreases with gt/a^ whereas the width 
AP decreases as 1/a. Some have gone as far as to define the weak value as 
^'^{jm') [11]- Nonetheless, we will show that (|P^ should not be inter- 
preted as an insignificant artifact of the weak measurement procedure and has 
an integral role in measuring the A'"th-order joint weak value. 

One can express the full weak value in terms of the two expectation values of 
the pointer. 



= -(-^ + ^tP) ■ (18) 
gt \2a h /fi ^ ^ 

In their derivation of weak values, AAV made the natural choice of a Gaus- 
sian for the initial pointer state, as do we. This state also happens to be 
the ground state |0) of a harmonic oscillator with mass m and frequency uj. 
For illustration, if one reparameterizes the width of the Gaussian in terms of 
muj such that a — ^Jh/2muj it becomes apparent that the operator in the 
expectation value in Eq. (18) is just the familiar lowering operator. 



—X + r.^P. (19) 

The operator in Eq. (18) will transform the pointer just as the lowering op- 
erator does, even though the pointer is not actually in a harmonic potential. 
This fact will simplify some of the following calculations. Furthermore, now 
the weak value can be re-expressed as: 

i\ =^(a),.. (20) 
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To our knowledge, this is the first time in the hterature that this simple but 
important relationship between the annihilation operator and weak measure- 
ment has been described. The reason the annihilation operator is related to 
the weak value can be understood as follows. When the couphng is sufficiently 
weak, the expansion in Eq. (6) shows that the largest pointer amplitude is left 
unchanged in the ground state. The interaction Hamiltonian shifts some of 
the pointer state into the first excited state by creating a small amplitude, 
proportional to gtA, for the |1) state. If we restrict ourselves to the post- 
selected subensemble, as in Eq. (9), then this small amphtude changes to be 
proportional to gt (^)^- The annihilation operator removes the part of the 
state that is left unchanged by the couphng, leaving only the shifted compo- 
nent. In other words, the annihilation operator isolates only that part of the 
pointer state that is changed by the interaction. 

We now move on to a derivation of A^th-order joint weak values. In this 
section, we combine the strategy outlined in the introduction for measuring 
joint observables with the use of the annihilation operator to extract the weak 
value. As in previous works, to measure the operator M = YijLi we apply 
separate von Neumann interactions coupling each Aj to its own pointer, 
as in Eq. (3) [22,23]. To simplify the expressions to come we set all Qj to 
be equal and rewrite the momentum operators Pj in terms of the respective 
raising and lowering operators, a] and a^ , for each of the pointers, 

n^i^j:A,{a]-a,). (21) 

Now we require N different pointers, all beginning in an initial state defined 
by Eq. (5). The total initial pointer state can be described by the ground 
state of N harmonic oscillators: 

N 

i^) = nK-) = ior''- (22) 

Continuing, using the number-state notation to describe the pointer, we cal- 
culate the state of the combined system after the interaction Hamiltonian is 
applied. 
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(23) 
(24) 

(25) 
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where is the state where the jth pointer is in the first-excited state and 
all the other pointers are in the ground state (e.g. IO1I2O3...OJV)). Here, we 
have expanded the state in powers of gt. Eq. (23) shows that to first order, 
the interaction Hamiltonian coupling the measuring device to the system can 
displace only one of the N pointers at a time. Simultaneous shifts of multiple 
pointers come from higher-order terms in the propagator. We are particularly 
interested in the A^th term in the expansion, 

This term is the lowest-order one in the expansion which can simultaneously 
transfer all N pointers into the first excited state (e.g. |lil2l3...lAr))- This 
state, which we label as , is created when each term in the above sum 

supplies one raising operator. The terms in the sum can contribute the 
distinct raising operators in any order and so the portion of Eq. (26) that 
creates the state is equal to 

where p denotes the sum of all A'"! orderings of the set of N operators 

|l/fc|. Note that these different orderings are only distinct when the operators 
do not commute. The remaining portions of Eq. (26) create states where at 
least one pointer is left in the initial state (e.g. |2i02l3...lAr)). Projecting onto 
{F\ completes the post-selection and leaves us with. 



{F\ oxp ( ^) \0f'- = lO)''-'^' {F\I) + g X: {F\ A,\I)\ I,) + ... (28 



We renormalize the resulting A'"-pointer state |$fi) and then truncate the 
amplitude of each term at the lowest nonzero order in gt, 

l$/.)-|o) ^ +••• 

(30) 

This is equivalent to dividing by {F\I), the renormalization constant in the 
limit of no coupling. In analogy with Eq. (20), we now wish to take the expec- 
tation value of the product of the annihilation operators for all N pointers. 
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N 

0=U &j. (31) 

In Eq. (30), the ll)*^^ state is the lowest order term that does not go to zero 
when acted on by O; this term becomes, 

Clearly, to lowest nonzero order the expectation value then becomes. 



'p)^={^f,\0\^f,) (33) 
1 {F\P{A,}JI)_^ ^^^^ 



2a) TV! 

The next lowest order term in the expectation value corresponds to any of the 

pointers undergoing an extra pair of transitions (i.e., a pointer is raised to 
|2) and subsequently lowered back to |1)). Consequently it will be reduced in 

size by a factor of 2 (|^) compared to the lowest order term. Using Eq. (35) 



,2cr ) 

the iVth-order joint weak value can now be expressed in the simple formula 




It is often the case that each operator Aj acts on a different particle, ensuring 
that all commute. This allows the further simplification of the A^th-order 
joint weak value to 




For commuting observables, the magnitude of the simultaneous shift in the 
N pointers that results from concurrent kicks from all N terms in the Hamil- 
tonian in Eq. (3) is proportional to the shift in one pointer created by a 
single von Neumann Hamiltonian for measuring operator M . The role of the 
annihilation operators is to isolate this simultaneous pointer shift from the 
total uncorrelated shifts of the N pointers and thus duplicate the action of 
7i = gMP, without the need for multiparticle interactions. 
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Since Eq. (37) requires the measurement of the annihilation operator, which is 
not Hermitian, one might think the expression is, in principle, unmeasurable. 
In fact, if one expands the annihilation operator in terms of X and P for each 
pointer then one is simply left with expectation values of products of X or P 
for each pointer. One then measures X in one ensemble of pointers and P in 
an identically-prepared ensemble. 

The expression in Eq. (37) is the central result of this work. As in previous 
papers, this result shows how one can practically measure a joint weak value 
even without the multiparticle interactions the AAV method requires [22,23]. 
However, this expression is much more elegant and it makes it clear that the 
annihilation operator plays a key role in joint weak measurements. Specifi- 
cally, with the use of the annihilation operator, the similarity to the strong 
measurement expectation value in Eq. (4) is apparent. For strong measure- 
ment, the equivalent expectation value to the A^th-order joint weak value is 




The similarity is striking and makes a good case for the use of the annihilation 
operator in the understanding of weak values. 

Lets compare Eq. (37) to the previous results for the A^th-order joint weak 
value [23]. In the previous paper, it was expressed recursively in terms of 
two A^th-order correlations between the pointers and to A^ different joint weak 
values of order A^ — 1. Utilizing this recursive formula, the A^th-order joint 
weak value can be expressed purely in terms of the expectation value of posi- 
tion and momentum correlations. This expression includes 2^"'"^ — 2 distinct 
correlations of various orders, although most will be close to the A^/2 order as 
the number of distinct expectation values at each order follows the binomial 
distribution. In comparison, Eq. (37) relates the A^th-order joint weak value 
to 2^ correlations in the positions and momenta of all A^ pointers and so re- 
quires roughly half the number of expectation values as the final result from 
the previous paper (but of higher order) . 

As a specific example of the use of Eq. (37), the weak value of the product of 
two spin components SixS2y would be. 



{SlxS2y)^-{^-^ {0,1(12) 



(39) 



'2aV 1(1'-^ \ / 1 .a 



The real and imaginary parts of the weak value are then 
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Re {S^S,,)„ = (i)' {{XA}^^ - ^ (AA) J (41) 
Im = (^)' ((^.ft)^, + (AX,)^J . (42) 

The importance of the pointer momentum shift is demonstrated in the above 
example. With our measurement technique even the real part of weak value is 
related to the pointers' momenta, Pi and P2 . In general, the momentum and 
position observables for each of the N pointers will appear in the expression 
for the real part of the A^th-order joint weak value. 

Note that like single weak measurements, this method for measuring the 
A^th-ordcr joint weak value is not limited to the particular interaction or 
pointer used in our measurement model [24]. For example, one can perform 
a derivation very similar to the one presented here where a spin, as opposed 
to position, pointer is used. For a spin pointer, the Hamiltonian would be 
H — —gASy — igA{S^ — S~)/2, where and S~ are the raising and lower- 
ing operators for the Si basis. The initial pointer state would be the lowest 
eigenstate of 5"^, with eigenvalue —hs. In this case, the expression for the 
A'^th-order joint weak value in terms of N spin pointers is 




where is the z-basis lowering operator for the jth pointer and all Aj are 
assumed to commute. An important advantage of using spin is the absence of 
unequal coefficients in the expression for the lowering operator. This puts the 
shifts in the pointer observable and its conjugate on equal footing. Using such 
a pointer means that the physical shift in the conjugate observable does not 
become smaller as the measurement becomes weaker. Expectation values are 
also particularly easy to measure for spins (and polarizations), especially spin 
1/2 systems since there are only two basis states which need to be projected 
onto. For instance, the iVth-order joint weak value would only require 2^^^ 
measurements in total if N spin 1/2 pointers were used. 

In the present work, we have greatly simplified a recent extension of weak mea- 
surement which makes the experimental investigation of composite, or joint 
observables possible [22,23]. We have shown that when single and joint weak 
values are expressed as expectation values of annihilation operators, they take 
on a surprisingly elegant form very similar to that seen in standard strong 
measurement. This form is easily generalized to any measurement device in 
which the initial pointer state is the eigenstate of an appropriate lowering op- 
erator. With the extension, the weak measurement of joint observables only 
requires the same apparatus that one would need to weakly measure each of 
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the component observables separately. Joint observables are central to the de- 
tection and utilization of entanglement in multiparticle systems. The weak 
measurement of these observables should be particularly useful for investigat- 
ing post-selected systems such as those that have been used produce novel 
multiparticle entangled states or those that implement quantum logic gates 
[25,26]. 
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